EPJ manuscript No. 

(will be inserted by the editor) 



Statistical Properties of Eigenvalues for an Operating Quantum 
Computer with Static Imperfections 

Giuliano Benenti 1 , Giulio Casati 1,2 , Simone Montangero 1 , and Dima L. Shepelyansky 3 



o 
o 

c 



> 
o 
m 

O 
(N 
O 

Oh 

-I— > 

C 

o3 
=3 
0< 



1 International Center for the Study of Dynamical Systems, Universita degli Studi dell'Insubria and Istituto Nazionale per la 
Fisica della Materia, Unita di Como, Via Valleggio 11, 22100 Como, Italy 

2 Istituto Nazionale di Fisica Nucleare, Sezione di Milano, Via Celoria 16, 20133 Milano, Italy 

3 Laboratoire de Physique Quantique, UMR 5626 du CNRS, Universite Paul Sabatier, 31062 Toulouse Cedex 4, France 

Received: June 19, 2002 

Abstract. We investigate the transition to quantum chaos, induced by static imperfections, for an operat- 
ing quantum computer that simulates efficiently a dynamical quantum system, the sawtooth map. For the 
different dynamical regimes of the map, we discuss the quantum chaos border induced by static imperfec- 
tions by analyzing the statistical properties of the quantum computer eigenvalues. For small imperfection 
strengths the level spacing statistics is close to the case of quasi-integrable systems while above the border 
it is described by the random matrix theory. We have found that the border drops exponentially with the 
number of qubits, both in the ergodic and quasi-integrable dynamical regimes of the map characterized by 
a complex phase space structure. On the contrary, the regime with integrable map dynamics remains more 
stable against static imperfections since in this case the border drops only algebraically with the number 
of qubits. 



PACS. 03.67.Lx Quantum computation 
Many-body theory 



05.45.Mt Semiclassical chaos ("quantum chaos") - 24.10.Cn 



1 Introduction 

Quantum computers, if constructed, would be capable of 
solving some computational problems much more efficiently 
than classical computers Jjj . Shor Q constructed a quan- 
tum algorithm which performs integer factorization expo- 
nentially faster than any known classical algorithm. It was 
also shown by Grover [0 that the search of an item in an 
unstructured list can be done with a square root speedup 
over any classical algorithm. These results motivated a 
great body of experimental proposals for the construction 
of a realistic quantum computer (see Q and references 
therein). 

While the technological challenge to develop scalable, 
fault tolerant quantum processors is highly demanding, it 
is well appreciated that decoherence, due to the coupling 
with the environment, will be the ultimate obstacle to the 
realization of such devices. In addition, even in the ideal 
case in which the quantum computer is isolated from the 
external world, a proper operability of the computer is not 
guaranteed: unavoidable internal static imperfections in 
the quantum computer hardware represent another source 
of errors. For example, the energy spacing between the 
two states of each qubit can fluctuate, e.g. due to mag- 
netic field inhomogeneities in nuclear magnetic resonance 
quantum processors jjj. Moreover, since qubit interactions 
are required to operate two-qubit gates and generate en- 



tangled states, unwanted residual interactions will appear. 
For example, when the inter-qubit coupling is switched 
off, e.g. via a potential barrier created by a point contact 
gate in the quantum dots proposal 0], some unavoidable 
residual interactions still remain. Therefore the quantum 
computer hardware can be modeled as a qubit lattice and 
one has to consider it as a quantum many-qubit (-body) 
interacting system Q . Many-body quantum systems have 
been widely investigated in the field of quantum chaos 
[HQ] and it is now well known that residual interactions 
can lead to quantum chaos characterized by ergodicity of 
the eigenstates and level s pacing statistics described by 
random matrix theory ]^j9|, ^0| , p]| . In the regime of quan- 
tum chaos the wave functions and energy spectra become 
so complicated that a statistical description should be ap- 
plied to them. Thus it is important to study the stability 
of quantum information processing in the presence of real- 
istic models of quantum computer hardware imperfections 
and in concrete examples of quantum algorithms. The sta- 
bility of the static quantum computer hardware was stud- 
ied in [^j and it was shown that residual inter-qubit inter- 
actions can destroy a generic state stored in the quantum 
computer if the amplitude of static imperfections is above 
the quantum chaos border. For a static quantum hard- 
ware this border drops linearly with the number of qubits 
while the average energy level spacing drops exponentially. 
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Above the quantum chaos border an exponential number 
of ideal multi-qubit states becomes mixed after a finite 
chaotic time scale. The implications of such imperfections 
for the Grover algorithm and quantum Fourier transform 
were analyzed in Ref. p2| . 

In this paper, we address the question of the transi- 
tion to chaos for an operating quantum computer, which 
is running an efficient quantum algorithm simulating the 
dynamical evolution of the so-called quantum sawtooth 
map The sawtooth map is a paradigm of classical 
]l4| and quantum Jl5| chaos, and exhibits a variety of 
different behaviors, from anomalous diffusion to quantum 
ergodicity and dynamical localization. We note that in the 
above quantum algorithm the number of redundant qubits 
can be reduced to zero and complex dynamics can be in- 
vestigated already with less than 10 qubits. It is there- 
fore important to study the stability of such algorithm, in 
view of its possible implementation in the first generation 
of quantum processors operating with a small number of 
qubits 0,0. 

Since spectral statistics is the most valuable tool to 
detect the transition to quantum chaos [||^] , we study the 
statistical properties of the eigenvalues of the quantum 
sawtooth map, simulating different dynamical regimes (er- 
godic, quasi- integrable, and integrable), in the presence of 
static imperfections in the quantum computer hardware. 
We will show that the threshold for the transition to chaos 
induced by static imperfections drops exponentially with 
the number of qubits both in the ergodic and in the quasi- 
integrable regime, while in the integrable regime the chaos 
border drops only polynomially with the number of qubits. 
We note that this paper complements Ref. |l^] , devoted to 
the study of the eigenvectors of the model in the ergodic 
regime. 

The paper is organized as follows. In Section 2 we dis- 
cuss the quantum algorithm simulating the quantum saw- 
tooth map model and its realization in the presence of 
hardware static imperfections. The effect of these imper- 
fections on the spectral statistics of the model is described 
in Sections 3,4, and 5, for the ergodic, quasi-integrable, 
and integrable regime, respectively. Our conclusions are 
presented in Section 6. 



2 The model and the quantum algorithm 

The quantum sawtooth map is the quantized version of 
the classical sawtooth map, which is given by 

n = n + k(9 -tt), 9 = 9 + Tn, (1) 

where (n, 9) are conjugated action-angle variables (0 < 
8 < 2tt), and the bars denote the variables after one 
map iteration. Introducing the rescaled momentum vari- 
able p = Tn, one can see that the classical dynamics de- 
pends only on the single parameter K — kT. As it is 
known, the classical motion is stable for —4 < K < and 
completely chaotic for K < — 4 and K > ]l4| . 

The quantum evolution for one map iteration is de- 
scribed by a unitary operator Uq (Floquet operator) act- 



ing on the wave function tp: 

1P = U^ = e -TftV2 e ,fe(e-) 2 /2 V)j (2) 

where h = —id/ 88 and ip(9 + 2tt) = ip(8) (we set Ti = 1). 
The classical limit corresponds to k — > oo, T — > 0, and 
K = kT = const. 

In this paper, we study the quantum sawtooth map (|^) 
closed on the torus — tt < p < tt. Therefore the classical 
limit is obtained by increasing the number of qubits n q — 
log 2 N (N is the total number of levels), with T = 2%/N 
(k = K/T, -N/2 < n < N/2). We consider the ergodic, 
quasi-integrable and integrable regimes. To this end, we 
take K = V2, K = —0.1, and K = —1, respectively. 

The quantum algorithm introduced in |fl3f simulates 
efficiently the quantum dynamics (|^) using a register of 
n q qubits. It is based on the forward/backward quantum 
Fourier transform ]l9| ] between the 9 and n representa- 
tions and has some elements of the quantum algorithm 
for kicked rotator J2(J. Such an approach is rather conve- 
nient since the Floquet operator Uq is the product of two 
operators U k = e ^e-^) 2 /2 and (j T = e -iTh 2 /2. the first 

one is diagonal in the 8 representation, the latter in the n 
representation. Thus the quantum algorithm for one map 
iteration requires the following steps: 

I. The unitary operator Uk is decomposed in n q two-qubit 
gates 

where 9 = 2TrJ2 a i^ 1 ^ with € {0,1}. Each two-qubit 
gate can be written in the {|00), |01), |10), |11)} basis as 
exp(ifc7r 2 Z?), where D is a diagonal matrix with elements 

{— --(--—) --(--—) 

X 2n\ n q \2i 2n J ' n q \2* 2n J ' 

\2 l 2n q J V2 J 2n q J f V ; 

II. The change from the 8 to the n representation is ob- 
tained by means of the quantum Fourier transform, which 
requires n q Hadamard gates and n q (n q — l)/2 controllcd- 
phase shift gates Jn]] . 

III. In the new representation the operator Ut has essen- 
tially the same form as Uk in step I and therefore it can 
be decomposed in n q gates similarly to equation ||). 

IV. We go back to the initial 9 representation via inverse 
quantum Fourier transform. 

On the whole the algorithm requires 3n q + n q gates per 
map iteration. Therefore it is exponentially efficient with 
respect to any known classical algorithm. Indeed the most 
efficient way to simulate the quantum dynamics (Q) on 
a classical computer is based on forward/backward fast 
Fourier transform and requires 0(n q 2 nq ) operations. We 
stress that this quantum algorithm does not need any ex- 
tra work space qubit. This is due to the fact that for the 
quantum sawtooth map the kick operator Uk has the same 
quadratic form as the free rotation operator Ut- 
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We model the quantum computer hardware as an lin- 
ear array of qubits (spin halves) with static imperfections, 
i.e. fluctuations in the individual qubit energies and resid- 
ual short-range inter-qubit couplings § . The model is de- 
scribed by the following many-body Hamiltonian: 



H 8 = J2(^o + S l )at 



2^ J ij°i°j 
i<3 



(5) 



where the &i are the Pauli matrices for the qubit i, and 
A is the average level spacing for one qubit. The second 
sum in (|^) runs over nearest-neighbor qubit pairs, zero 
boundary conditions are applied, and 8i, Jij are randomly 
and uniformly distributed in the intervals [—8/2, 8/2] and 
[— J, J], respectively. We model the implementation of the 
above described algorithm as a sequence of instantaneous 
and perfect one- and two-qubit gates, separated by a time 
interval r g , during which the hardware Hamiltonian (^J) 
gives unwanted phase rotations and qubit couplings. There- 
fore we study numerically the many-body Hamiltonian 



H(T)=H s + H g ( T ), 



where 



kr g )h k . 



(6) 



(7) 



Q_ 




Fig. 1. The level spacing distribution for K — y/2, n q = 11 
and J = 0, at e = (circles) and e = 10 -3 (squares). The 
dashed and full lines give the theoretical distribution (hot) and 
([ill), respectively. Inset: statistics for one symmetry class at 
e = and same parameter values as in the main figure. The 
dashed line is the theoretical distribution (^). 



Here hk realizes the fc-th elementary gate according to the 
sequence prescribed by the above algorithm. We assume 
that the average phase accumulation given by An is elim- 
inated, e.g. by means of refocusing techniques fl2l|| . 

Since the evolution operator (|2j) remains periodic in 
the presence of static imperfections, C/ e . p (r + T) = U^ p (r) 
(e = 8t 9 , p = Jr g rescaled imperfection strengths), all 
information about the system dynamics is included in the 
quasi-energy eigenvalues Xa ' P ' and eigenstates <f>a of the 
perturbed Floquet operator: 



UeATm 



exp 



(iX^)&»\ 



(8) 



We study numerically the spectral statistics of the Flo- 
quet operator (^) for a quantum computer running the 
quantum sawtooth map algorithm in the presence of static 
imperfections described by (^). We construct numerically 
the Floquet operator in the momentum representation (i.e. 
the quantum register states basis) using the fact that the 
one step map evolution (including static imperfections) of 
each quantum register state gives a column in the matrix 
representation of this operator. 

We consider 4 < n q < 12 qubits. In order to reduce 
statistical fluctuations, data are averaged over 3 < No < 
10 3 random realizations of static imperfections. In this 
way the total number of Floquet quasi-energies is Nd x 
2 n " sa 10 4 , which is large enough to get stable results in 
the study of spectral statistics. 



3 The ergodic regime 

We focus our attention first on the ergodic regime, where 
the eigenfunctions of the unperturbed Floquet operator 



are given by a complex superposition of order N = 2 Tlq 
quantum register states. In particular, we consider K = 
\[2, where the corresponding classical sawtooth map ([j]) 
is completely chaotic. In the ideal case in the absence of 
static imperfections (e = 0, p = 0), the eigenvalues of the 
Floquet operator are divided in two different symmetry 
classes, even and odd with respect to the transformation 
n — > —n. Thus the Floquet operator has two blocks which 
can be studied independently. A convenient tool to char- 
acterize the spectral properties of the system is the level 
spacing statistics P{s), giving the probability to find two 
consecutive eigenvalues whose energy difference, normal- 
ized to the average level spacing, is in [s, s + ds]. Since we 
are in the ergodic regime, the level spacing statistics for 
e_ach block is well described by the random matrix theory 
in the presence of time-reversal symmetry: 



Po(s) 



(9) 



This theoretical distribution is in agreement with the nu- 
merical results shown in the inset of Fig. [I] for the symmet- 
ric class. The global spectral statistics can be computed 
from the single spectral statistics of each symmetry class 
||, and is given by: 



4/4 



( 10 ) 

again in good agreement with the numerical data of Fig. 

Since static imperfections break the time reversal sym- 
metry, the two symmetry classes become mixed and the 
system undergoes a crossover from (|Io| ) to a different ran- 
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dom matrix statistics 



32s 2 



-4s 2 /tt 



(11) 



The good agreement between this limiting distribution 
and the actual statistics computed from numerical data 
is shown in Fig. [l] for n q — 1 1 and e = 10~ 3 , J — 0. A con- 
venient quantity to characterize the crossover from one 
limiting distribution to another is pfl 



j;°{p{ s )-p u { s ))ds 
s; a {p ( o\ S )-Pu{ S ))ds' 



(12) 



where sq — 0.50285... is the first intersection point of 

(2) 

Pq (s) and Pu(s)- This parameter goes from one to zero, 
when the spacing probability changes from ([n]) to (O). 
The behavior of 77 as a function of e is shown in Fig7^|, 
for different n q values. This figure shows that it drops to 
zero when the imperfection strength grows. For e — > and 
small n qi there are significant deviations from the value 
77 = 1 corresponding to the limiting distribution ( |Io|) . This 
is due to the fact that the e = statistics includes only 
2 n q levels spacings. For very small e values the statistics 
remains poor since each static imperfections realization 
gives essentially the same P(s) distribution. Nevertheless, 
it is clear that r](e = 0) goes to one when the number of 
qubits increases. 

In order to study the dependence of the 77-crossover on 
the number of qubits, we define the critical imperfections 
strength e x as the e value at which r\ = 0.2 (similar results 
are obtained for different ij values) . In Fig. we show that 
e x drops exponentially with the number of qubits. This ex- 
ponential threshold is due to the fact that, in the ergodic 
regime, quantum eigenstates are given by a complex su- 
perposition of an exponentially large number of quantum 
register states. Indeed, due to quantum chaos, the eigen- 
states of the unperturbed (e = 0, J=0) Floquet operator 
(§|) can be written as 



N 



L 0) = E c « 



m) 



(13) 



Here u rn are the quantum register states, and are ran- 
domly fluctuating components, with \ca^\ ~ 1/V~N due 
to wave function normalization. The transition matrix ele- 
ments between unperturbed eigenstates can be computed 
in perturbation theory. For J — 0, they have a typical 
value 



^yp~|(^ 0) |E^^I^ 0) )l 
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Fig. 2. Dependence of the r\ parameter on the scale imper- 
fection strength e at K = v2, J = 0. From right to left: 
n q = 4,5, 12. 



In this expression, the typical phase error is S^/n^ (sum 
of n q random detunings <5j's) and T g n g ~ T g n q is the time 
used by the quantum computer to simulate one map step 
The last estimate in (n4|) results from the sum of ./V terms 



of amplitude \c a 



( m )J,f 



1/N and random phases. Since 



the spacing between quasi-energy eigenstates is AE ~ 
1/N, the threshold for the breaking of perturbation theory 
can be estimated as 



V typ /AE 
The analytical result 



e x n 5 q / 2 VN 



n\' 2 VN 



(15) 



(16) 



is confirmed by the numerical data of Fig. [| For the case 
J = 5, the threshold e x approximately decreases by a 
factor 1.2 with respect to the J = case (see again Fig. ||), 
since residual inter-qubits interactions introduce further 
couplings between Floquet eigenstates. However, the same 
functional dependence ( |l6|) takes place. 

We note that the same exponential sensitivity to static 
imperfections was detected for the Floquet eigenstates in 
Ref . jl8| . This confirms that spectral statistics are a useful 
tool to characterize the mixing of unperturbed eigenstates 
in an operating quantum computer. 



4 The quasi-integrable regime 



»=i 



i\ E 



c (m) c (n) 



m,n— 1 i— 1 



en" 



f 1 e <£ m) 4 m) *i ~ *< ,2 N- l/2 - 



m—l 



We now focus our study on the regime of quasi-integrability 
of the sawtooth map, in which the system is not chaotic 
* di(u n \a*\u m )\ (14) (i.e. the Lyapunov exponent is zero) and there is a non in- 

tegrable component in the phase space. In this situation, 
the system exhibits interesting physical properties, such 
as anomalous diffusion and hierarchical phase space struc- 
tures IlS]. In this Section, we consider the case K = —0.1 
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Fig. 3. Dependence of the imperfections strength e x at which 
r\ = 0.2 on the number of qubits, for K = v2, J = (di- 
amonds) and J = 8 (circles). The line gives the theoreti- 



cal dependence e x = A2 
A = 4.3. 



3 , with the fitting constant 



and we study the transition to chaos for the eigenstates 
of a quantum computer simulating the quantum sawtooth 
map. 

A useful tool to demonstrate the transition to chaos 
induced by static imperfections is the parametric depen- 
dence of the quasi-energy eigenvalues. The evolution of 
a part of the spectrum as a function of the imperfection 
strength e (at J = 0) is shown in Fig. ||. It makes evident 
the qualitative change induced by static imperfections. At 
small e the spectrum exhibits quasi-degeneracyes, while 
at large e avoided crossings appear. Indeed, at small e, 
Floquet eigenstates with very close eigenvalues may lay 
so far apart that their overlap is negligible. Thus there 
is essentially no level repulsion for these eigenvalues. On 
the contrary, at large e, Floquet eigenstates are delocal- 
ized and therefore their overlap is significant, and induces 
level repulsion |22|. The delocalizing effect of static im- 
perfections is evident in the Husimi functions of Fig. [5], 
drawn from two typical Floquet eigenstates |23[|24|| . 

The two top plots of Fig. || represent two exact eigen- 
states (e = 0, J = 0). In the right picture, the quan- 
tum probability is concentrated around an ellipse corre- 
sponding to a classical integrable trajectory (a torus in the 
phase space (0,p)). We note that the map (||) is the dis- 
cretized time evolution for an harmonic oscillator. There- 
fore it gives elliptic trajectories as far as border effects 
can be neglected. A completely different kind of eigenvec- 
tor appears on the left picture. The Husimi function is 
spread, and this reflects the properties of the correspond- 
ing classical non-integrable trajectories which diffuse in a 
non-Brownian way |13| . We note that the number of such 
eigenstates is non- negligible. Indeed, we have checked that 
the probability of finding a diffusive classical trajectory 
starting from a randomly chosen initial condition is 0.12. 
We also note that classical diffusion is suppressed by quan- 
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Fig. 4. Dependence of quasi-energy eigenvalues on the imper- 
fections strength e, for n q = 9, K = —0.1, J = 0. The thick 
dashed (dot-dashed) curves corresponds to the eigenvector rep- 
resented in Fig. ^ right (left). 



turn interference effects giving a distribution localized in 
momentum. The remaining pictures of Fig. [| represent 
the same Floquet eigenstates, computed in the presence 
of static imperfections. The middle figures are obtained 
slightly below the chaos border induced by static imper- 
fections, and maintain same similarities with the exact 
eigenstates. On the contrary, in the bottom pictures, taken 
above the chaos border, the eigenfunctions are spread in 
all the available phase space and any structure has been 
destroyed. We note that, starting from completely differ- 
ent unperturbed eigenfunctions, one gets statistically in- 
distinguishable eigenfunctions with randomly fluctuating 
components. 

We now characterize the transition by studying the 
spectral statistics. In this regime, the limiting P(s) dis- 
tribution at e = is non trivial. Indeed, the Shnirelman 
theorem states that for a nearly integrable system the level 
spacing statistics P(s) exhibits a huge peak near the ori- 
gin (s = 0) p5| , p6[ . In the following we choose to eliminate 
such peak, related to time-reversal invariance [n — > —n, 
9 — ► 2ir — 0). To do that we operate the transformation 



(17) 



where <j> plays the role of an Aharonov-Bohm flux p8[ . 
The spectral statistics is still complex since in the classi- 
cal limit the phase space has two components, integrable 
and non integrable. In the semi-classical limit, the levels 
belonging to these two components become uncorrelated 
and the global spacing statistics is given by the superpo- 
sition of each component's statistics |^9|. We note that in 
this case both components have non-trivial statistics. In 
particular it has been shown that for the related harmonic 
oscillator case the quasi-energy spectral statistics exhibits 
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Fig. 5. Husimi functions in action-angle variable (p, 6) (—tv < 
p < 7r -vertical axis- and < 6 < 2tt -horizontal axis-) for the 
quantum sawtooth map at n q = 9, K = —0.1, J — 0. The left 
(right) column corresponds to the dot-dashed (dashed) line in 
Fig. [§ for £ = (top), e = 10" 3 (middle), e = 3 x 10" 3 
(bottom). We choose the ratio of the action-angle uncertainties 
s = Ap/A9 = TAn/AG = 1 (A P A6 = T/2). The color is 
proportional to the density: blue for zero and red for maximal 
density. 

level repulsion and a finite number of peaks j3Cj, instead 
of the Poisson statistics typical of integrable systems. 

The spectral statistics at e = 0, J = is shown in Fig. 
|(| for n q = 11. In the same figure, one can see that static 
imperfections induce a crossover to the Wigner-Dyson dis- 
tribution (|TT|). One can compare this figure with Fig. [p. 
The e = distributions are completely different, reflect- 
ing different dynamical regimes. On the contrary, the same 
universal distribution is found in the regime in which static 
imperfections destroy all symmetries. Even though in this 
case we cannot provide an analytical expression for the 
e = J = statistics, the crossover to the Wigner-Dyson 
distribution (jll]) can be characterized by the parameter 

fj=U [P( S )-Pu( S )} 2 dsj , (18) 

measuring the distance of the spacing distribution from 
(|ll|). The behavior of fj as a function of e (at J = 0) is 




Fig. 6. Level spacing distribution for n q = 11, K = —0.1, 
8o = 4> = v2/5, J = 0, e = (circles) and e = 10~ 3 (squares). 
The solid line gives the Wigner-Dyson distribution (Eq. (|ll|)). 




8 



Fig. 7. Dependence of the parameter fj on the imperfection 
strength e for K = -0.1, 8 = <t> = J = 0. From right to 

left: n g = 4, 5, . . . , 12. 



shown in Fig. |7] for various n q values. This figure shows 
again that the Wigner -Dyson distribution (fj = 0) is 
reached by increasing e. 

Similarly to what we have done in the ergodic regime, 
we define a critical imperfections strength e x as the e value 
at which fj — 0.2. The dependence of e x on the number of 
qubits is shown in Fig. ||. Quite surprisingly, even in this 
quasi-integrable regime we can fit the exponential drop of 
e x with n q via the analytical function e x = B2~ nq / 2 n q 5 ^ 2 , 
with the fitting constant B = 6.5. Therefore, it seems that 
the theoretical argument developed in the previous Sec- 
tion for ergodic eigenstates has a broader validity and can 
apply also in the more general mixed phase space dynam- 
ics. 
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Fig. 8. Dependence of the imperfections strength e x at which Fig. 9. Quasi-energy levels X a at n q = 9, K — — 1, e = J = 



r\ = 0.2 (see Fig. (?]) on the number of qubits n q (diamonds), 
for K — —0.1, do = 4> = v2/5, J = 0. Squares give the 
critical imperfection strengths at which the quantum eigen- 
states entropy 5 = 1. Lines give the theoretical dependency 
e x = BN' 1/2 nq 5/2 , with the fitting constant B = 1 (below) 
and B — 6.5 (above). 



The same exponential sensitivity to static imperfec- 
tions can be detected also in the Floquet eigenstates. The 
mixing of unperturbed eigenstates, induced by static im- 
perfections, is characterized by the quantum eigenstates 
entropy 

N 

S a = -y^p Q/ 3l0g 2 p Q/3 , (19) 



(thick bands) and theoretical values (dashed lines). 
10" 1 



W 10" 2 



0=1 



where p a p — |(<^|</>a' p ' ) )| 2 . In this way S a = if 



0,p) 

a 

coincides with one eigenstate at e = p = 0, S a = 1 if <j>a ,p ^ 
is equally composed of two ideal (e = p — 0) eigenstates, 

and = n 9 (maximal value) if all <j)^ (/3 = 1, ...,7V 



10 




2 n «) contribute equally to 0q' p ' ) . In Fig. || we show that 
the critical imperfection strength at which S = 1 drops 
exponentially with the number of qubits (S is the average 
of S a over a and static imperfection realizations). This 
shows that the transition to chaos can be characterized 
both by the mixing of unperturbed eigenstates and by the 
transition in the spectral statistics. 



5 The integrable regime 

We now consider the integrable case K = — 1. Indepen- 
dently of n q , the quasi-energy spectrum is composed of 

6 degenerate levels (see Fig. 9), = (2ir/&)j (up to an 
unessential global phase factor). This phenomenon can be 
explained from classical mechanics: indeed the 6-th iterate 
of the map (0) gives the identity, i.e. all orbits are periodic 
with period at most 6 f3l| . The same phenomenon can be 
observed at K — — 2 (K = —3), where the spectrum has a 



Fig. 

fj = 



10. Dependence of the imperfections strength e x at which 



0.2 on the number of qubits n q , for K 



-1, J 







(diamonds) and J = S (circles); note the log-log scale. Lines 



give the theoretical dependence e x 



Cn, 



-5/2 



with the fitting 



constant C = 1 (below) and C — 1.8 (above). 



degeneracy 4 (3) and the 4-th (third) iterate of ([!]) is the 
identity. In these cases, the quantum evolution and the 
classical discretized evolution coincide: this means that 
the quantum unitary evolution (^|) maps a given phase 
space distribution in the same way as the discretized Li- 
ouville operator does in classical mechanics (see, e.g., Ref. 

§!)■ 

Static imperfections induce a crossover to the universal 
Wigner-Dyson distribution (|Tl"| ) , that can be characterized 
by the parameter fj defined in Eq. (118). Again, one can 



define a critical imperfection strength^ e x as the value of 
e at which fj(e x ) — 0.2. The result is shown in Fig. 10, 
where it is seen that e x drops polynomially with n q : e x oc 

n q . This happens both for J = and J — S. We note 
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that this algebraic dependence contrasts the exponential 
decay observed in previous Sections for the ergodic and 
quasi-integrable regimes. This can be understood via the 
following argument. Quantum ergodicity can be reached 
only when levels of the different bands of Fig. 9 are mixed. 
Non-degenerate levels are separated by an energy spacing 
AE which is ^-independent. On the other hand, also the 
typical overlap between nearby eigenstates (in space) is 
iV-independent (their space separation drops with N, but 
also their typical width drops with N). Since the typical 
error is 5*Jn^ and the time needed to simulate one map 
step is ~ T g riq (see the discussion following Eq. (|l4])), one 
can estimate the typical transition matrix element V typ ~ 

5/2 

en q , and the threshold for the breaking of perturbation 
theory as 



V typ /AE ~ e x n 5 J 2 



giving the analytical estimate 



-5/2 



(20) 



(21) 



6 Conclusions 

In this paper, we have studied the transition to quan- 
tum chaos, induced by static imperfections, in the spectral 
statistics of an operating quantum computer. The thresh- 
old for the transition to chaos drops exponentially with 
the number of qubits, both in the ergodic and in the more 
general quasi-integrable regime. On the contrary, in the in- 
tegrable regime the chaos border drops only algebraically 
with the number of qubits. This is due to the presence 
of a finite number of bands in the spectrum of the time 
evolution operator related to global periodicity of classi- 
cal dynamics. We note that similar strong spectral de- 
generacyes have been observed in the Grover's algorithm 
and in the quantum Fourier transform in Rcf. p2| . We 
also stress that a complex dynamical system is generically 
quasi-integrable. The simulation of such systems is well 
accessible to the first generation of quantum computers 
with less than 10 qubits and we think that this class of 
quantum algorithms deserves further studies. 
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